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^ ■ A study of saturation effects in two-hard-scale hadronic processes such as Mueller- 

Navelet jets is presented. The cross-sections are expressed in the dipole framework 
while saturation is implemented via an extention of the Golec-Biernat and Wiisthoff 
model. The transition to saturation is found to be more abrupt than in 7*— 7* 
cross-sections. Observablcs with a potentially clear saturation signal are proposed. 

1 Introduction 

^O i Hard processes involving two perturbative scales lead to cross-sections whose lin- 

^P ■ ear high-energy behavior is described by the well-known Balitsky-Fadin-Kuraev- 

^T I Lipatov (BFKL) PP equation. However, to respect whatever constraints unitarity 

may impose, it is well-known that the BFKL equation has to be modified beyond 

some energy limit, in order to describe cross-sections that saturate. Physically, the 

idea is that the gluon density in the BFKL ladder grows higher as one increases 

^ il the energy and that eventually recombinations will occur, limiting the number of 

^ • gluons in the ladder. 

Three measurements for studying this behavior can be considered: the 7*^7* 

total cross-section in e+e~ scattering |2I3| . Mueller-Navelet jets in hadron-hadron 

collisions ^, and forward jets in deep inelastic scattering |5I6) . The perturbative 

«2j ' scales in these processes are set by either the virtualities of the reaction-initiating 

C^ I photons or the transverse momenta of the measured jets. The aim of this work is 

to describe in a simple way how saturation effects could appear in those processes. 
Following the approach of Golec-Biernat and Wiisthoff whose saturation model 
|7] (GBW) for the proton structure functions provides a simple and elegant for- 
mulation of the transition to saturation, we will implement saturation effects in 
the dipole framework |8l9j . The basis of this approach is to consider that the inci- 
dent particules fluctuate into colorless quark-antiquark pairs (dipoles) which then 
interact. Saturation will then be modeled through the dipole-dipole scattering. 

While such a study has already been done for the 7*— 7* cross-section ^U], our 
work will focus on Mueller-Navelet jets; the extension to the forward-jet case will 
be straightforward. The key difference between the 7*— 7* and the Mueller-Navelet 
jet measurements is that the hard probes are no more virtual photons but the final- 
state jets. The functions expressing the fluctuation of a virtual photon into dipoles 
are know from QED, but the descrition of a forward jet in terms of dipoles requires 
more care. A first part is devoted to this problem and then saturation predictions 
within the GBW model are presented. Observables to be studied are proposed. 
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2 Formulation 

Mueller-Navelet jets are processes in which a proton strongly interacts with another 
proton or antiproton and where a jet with transverse niomemtum larger than a 
perturbative scale is detected in each of the two forward directions. Such hard 
processes obey the collinear factorization which allows one to deal only with hard 
cross-sections. The two cuts on the jets transverse momenta will be denoted Qi and 
Q2 and taken of the same magnitude in order to suppress the DGLAP evolution 
in the gluon ladder. The rapidity interval between the two jets Ary is taken to be 
large in order to lie in the high-energy regime. 

Considering first the leading-logarithmic approximation when the evolution is 
linear, the dipole formulation of this hard total cross-section reads: 

o-(Qi,Q2, A77) = / (fridzi (fr2dz2 (f>j{ri, zi,Ql) 0,/(r2, Z2, Q2) o-^d (^1 , ^2 , Ar/) , 

(1) 
where ri=i^2 are the transverse sizes of the dipoles interacting and 2:^=1.2 are the 
fractions of longitudinal momentum of the quarks in each dipole. The BFKL dipole- 
dipole cross-section is 

2i7r 7^(1 — 7) 

(2) 
where ■0(7) is the logarithmic derivative of the Gamma function. The dipole dis- 
tributions describing the forward-jet emissions have been denoted (f>j{ri, Zi, Qf)- 

Let us recall how one can obtain this dipole distribution. The fc^-factorization 
property JJ provides the general formalism for coupling external sources to the 
BFKL kernel through the convolution of impact factors. It can be proved that 
fcy-factorization is equivalent jl2| to the dipole factorization expressed by formula 
(^. The dipole distribution 4>j can thus be derived from the corresponding impact 
factors: the derivation [13114115] is made using the example of a final-state gluon 
with transverse momentum larger than Q being emitted off a perturbative onium 
{qq state) of size rg ^ 1/Kqcd- QCD factorization will allow to extend the result 
to the case of an incident hadron since the onium structure function factorizes out. 
Using fcr-factorization in the BFKL framework, the impact factor /(fc^,ro) of the 
onium-|-jet system is related to the elementary gluon-dipole coupling f'^(k'^,r) in 
the following way |14I15| : 

/(fc2,ro) = l^^f^log J-logQroj f d\ -^ MQr) f{k\r) (3) 

L TT xj ) J 2-nr 

in the collinear approximation Qtq 3> 1 for the onium. k is the transverse mo- 
mentum of the gluon connected to the BFKL kernel and x,j is the fraction of 
longitudinal momentum of the jet with respect to the onium. Formula Q can 
be interpreted as the equivalence for forward jets between the momentum-space 
(partonic) and coordinate-space (dipole) representations. The factor in brackets 
{(2asA''c/7r)logQrologl/xj} corresponds to the probability of finding a dipole 
of size 1/Q inside the onium of size tq; thanks to QCD factorization properties, 
it is included in the gluon structure function of the incident particule (here the 
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onium). f'^{k^,r) = (1 — Jo(fcr))/A:^ is nothing else than the gluon density inside 
the dipole of size r and, in the dipole formulation Q, is included in the dipole- 
dipole cross-section a^J . Having factorized out both the contribution to the struc- 
ture function and to the dipole-dipole cross-section, one is left with the function 
ff'ji'fjQ^) = J dz (j)j{r,z,Q^) which describes the resulting size distribution of the 
interacting dipole. Hence, one is led to identify 

^j{r,Q'') = :^Ji{Qr). (4) 

Zirr 

Let us now consider saturation effects. Initially, the GBW approach I?- is a 
model for the dipole-proton cross-section which includes the saturation damping 
of large-dipole configurations. For the description of 7*— 7* cross-sections at LEP 
|10|. it has been extended to dipole-dipole cross-sections: 

a,,(ri,r2,Ar;)=ao|l-exp(^-^^^)| . (5) 

The dipole-dipole effective radius r'^g{ri,r2) is defined through the two-gluon ex- 
change: 

2nalrl^{r„r,) ^ 4°i(ri, r„ 0) ^ 2nal min(r^ r,^) |l + log "^^^^'^^^^^ \ (6) 

while for the saturation radius i?o(A?7) = e^^i^'^-^va) /Q^ -^^6 shall use the same 
set of parameters as those in |7I1U| . that is A = 0.288, A770 = 8.1 for Qo = 1 GeV. 
Two other scenarios for reff (ri, r2) have also been considered: ^'cg = ^i''2/(''i + '^i) 
and r^fj = min(r^, r|). 

We shall use add in the hard cross-section ^ instead of a^^J to implement 
saturation in a simple way. However, in order to do so, one makes the non-trivial 
assumption that the dipole factorization still holds when the dipole-dipole cross- 
section is modified by saturation. 

3 Phenomenology 

Inserting Q and |SJ| in formula ^ leads to the simple final result for the Mueller- 
Navelet hard cross-sections modified by saturation within the GBW model: 

a(g„Q2, Ary)/ao = 1~2RIQ,Q2 / du'- -^—- /, ^^f^ 

(7) 
Some comments are in order. The dipole distribution 0,/(r, Q^) is not everywhere 
positive and we interpret this feature as a breakdown of the collinear approximation. 
It also means that one has to check that replacing a^J by add in O does not alter 
the positivity of the hard cross-sections, and this is indeed the case. Another check 
that our approximations require is that the cross-sections add '^ (^ofgff / ^^Rq^At]) , 
corresponding to the limit of small dipole sizes in |(SJ), lead to hard cross-sections be- 
having like 1/ {RqIAtj) max {Ql, Q2)} , as expected from transparency. The model 
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Tgfj = min(r^,r2) does not and therefore we cannot consider it in our approxima- 
tions. 

Let us investigate the phenomenofogical outcome, for hadron colhders, of our 
extension of the GBW models to Mueller-Navelet jets. The theoretical hard cross- 
sections are obtained from formula Q in terms of the physical variables Qi, Q2 
and A77. When plotting them, one observes the expected trend of the GBW model, 
that is a convergence of the cross-sections towards the full saturation limit cr -^ ctq . 
In order to appreciate more quantitatively the influence of saturation, it is most 
convenient to consider the quantities TZi/j defined as 

^ o-(Qi,Q2,A?7,) 

i. e. the cross-section ratios for two different values of the rapidity interval. These 
ratios display in a clear way the saturation effects. They also correspond to possible 
experimental observables since they can be obtained from measurements at fixed 
values of the jets longitudinal momenta xj^ and xj^ and thus are independent of the 
structure functions f{xj.,Q^) of the incident hadrons. Indeed, the experimental 
measurement is 

-^^^^ ^ fixj,,Ql)f{xj„QJ) aiQuQ2,A7^) (9) 

and the ratio of these cross-sections gives access to TZ. Such observables have 
actually been used for a study of Mueller-Navelet jets for testing BFKL predictions 
at the Tevatron |16I17| . 

In Fig.l we plot the values of TZ4.g/2A (resp. 72-8/4) as a function of Qi =Q2 = Q. 
72-4.6/2.4 is the observable that has been considered for the Tevatron |1()I17J while 
72-8/4 corresponds to realistic rapidity intervals for the LHC. As expected from the 
larger rapidity range, the decrease of 72 between the transparency regime and the 
saturated one is larger for the LHC than for the Tevatron. The striking feature of 
Fig.l is that the effect of saturation appears as a sharp transition for some critical 
range Q^1/Rq. No saturation effects would correspond to 72 constant equal to the 
high— Q^ limit of the plots while the full-saturated limit is 72=1. Comparing these 
ratios for Mueller-Navelet jets to those for the 7*— 7* measurement for the same 
values of the rapidity ranges, one interestingly sees that the 7*— 7* transition curve 
is much smoother, a phenomenon explained by the different structure of the dipole 
distributions. Indeed the formula to compute the 7*— 7* case is also formula Q) 
but with of course the well-known photon dipole distributions (j)^ instead of 4>j. As 
discussed in ^S], the dipole distribution 4>j{r, Q^) has a tail extending towards large 
dipole sizes, which are more damped by the saturation corrections. Hence (pj is 
more abruptly cut by saturation than the photon dipole distribution 0,^ . Note that 
saturation effects in forward-jets |18| can be studied in a straightforward manner 
using our formalism: it requires to combine Q with both dipole distributions (/)-y 
and 4)j. 

The signal displayed in Fig.l shows a clear transition to saturation, however 
the values of Q at which it occurs are rather low, probably to low for experi- 
mental Et — cuIs on jets. An interesting way out of this problem could be that 
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R8/4 and R4.6 /2.4 
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Figure 1. Cross-section ratios 'R-i/j- The resulting ratios for the two-gluon exchange model (1) 
and for r^g = rfr2/{ri +^2) (2) are plotted for rapidity intervals i = 8, j = 4 (highest set of curves) 



and i = 4.6, j = 2.4 (lowest set of curves). The comparison is made with 7* 
and equivalent kinematics. 



ratios for model 2 



the saturation scale is higher than the one we used in the present work, namely 
the one extracted from F2. Indeed, it has been proposed JR that the saturation 
scale could be higher for two-hard-scale processes like Mueller-Navelet jets than 
for one-hard-scale measurements like the proton structure functions. That would 
shift the transition shown in Fig.l towards higher Q. Another alternative to solve 
this "low-Q" problem would be to consider the detection of heavy vector or heavy- 
flavored mesons as alternatives to forward jets. Indeed, using J/'^'s, T's, D*'s, 
or B— mesons may provide hard probes of lower transverse momenta than jets, 
allowing to look deeper in the saturation regime. 

These possibilities of realizing hard hadronic probes of saturation certainly de- 
serve more studies in the near future. On the theoretical side, going beyong our 
approximations seems necessary while on the phenomenological side, simulations 
at Tevatron and LHC energies will be needed to give a quantitative estimate of the 
potential of hadronic colliders to reveal those new features of saturation. 
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